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Abstract 

Background: Identifying individuals at increased risk for melanoma could potentially improve public health through 
targeted surveillance and early detection. Studies have separately demonstrated significant associations between 
melanoma risk, melanocortin receptor (IVICIR) polymorphisms, and indoor ultraviolet light (UV) exposure. Existing 
melanoma risk prediction models do not include these factors; therefore, we investigated their potential to improve the 
performance of a risk model. 

Methods: Using 875 melanoma cases and 765 controls from the population-based Minnesota Skin Health Study we 
compared the predictive ability of a clinical melanoma risk model (Model A) to an enhanced model (Model F) using receiver 
operating characteristic (ROC) curves. Model A used self-reported conventional risk factors including mole phenotype 
categorized as "none", "few", "some" or "many" moles. Model F added MCIR genotype and measures of indoor and 
outdoor UV exposure to Model A. We also assessed the predictive ability of these models in subgroups stratified by mole 
phenotype (e.g. nevus-resistant ("none" and "few" moles) and nevus-prone ("some" and "many" moles)). 

Results: Model A (the reference model) yielded an area under the ROC curve (AUC) of 0.72 (95% CI = 0.69, 0.74). Model F was 
improved with an AUC = 0.74 (95% CI = 0.71-0.76, p<0.01). We also observed substantial variations in the AUCs of Models A 
& F when examined in the nevus-prone and nevus-resistant subgroups. 

Conc/usions: These results demonstrate that adding genotypic information and environmental exposure data can increase 
the predictive ability of a clinical melanoma risk model, especially among nevus-prone individuals. 
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introduction 

With a dismal overall survival for advanced disease, cutaneous 
melanoma results in more years-of-life-lost than any major adult 
cancer besides breast [1]. Fortunately, melanoma is highly curable 
when diagnosed at its earliest stages [2] , and a recent population- 
based skin cancer screening program in Germany successfully 
demonstrated a decrease in melanoma mortality [3] ; however, due 
to the relative rarity of melanoma compared to other cancers, like 
breast and prostate, implementing a population-based screening 
initiative can be expensive and impractical. Thus, refining such 
programs to prioritize patients at highest risk for the disease may 
improve the cost-effectiveness of such activities [4] . 



Although well-established melanoma risk factors exist (e.g. fair 
pigmentation, increased numbers of nevi, excessive outdoor 
ultraviolet (UV) light exposure), analysis of case-control studies 
demonstrate that as many as 50% of melanoma patients lack these 
risk factors, hindering our ability to identify those at greatest risk 
[5-7]. Additionally, certain relevant risk factors may be more 
common in different patients who develop different melanoma 
subtypes (i.e. superficial spreading, nodular, lentigo maligna, etc.). 
For example, patients who develop superficial spreading melano- 
ma are more likely to have an increased number of nevi, while 
lentigo maligna melanoma is more common in individuals with 
chronic sun damage. These clinical patterns have been synthesized 
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into a unifying concept suggesting that melanoma may result from 
difierent causal pathways (i.e. The Divergent Pathway Theory). 
Specifically, these studies describe two causal pathways of 
melanomagenesis associated with different UV exposure patterns 
and patient nevus phenotypes [8], the latter of which are known to 
have a strong genetic component [9,10]. This may explain why 
existing melanoma risk assessment tools have failed to reach the 
mainstream clinical arena, and suggests the need to account for 
this disease heterogeneity to maximize their accuracy. 

Fortunately, new risk factors have been described, including 
genetic markers and indoor tanning exposure, which may improve 
early detection and screening efforts [11-15]. Indoor UV 
exposure, which is associated with a 59% increase in melanoma 
risk when exposure occurs before the age of 35 [16,17], is of 
particular concern as approximately 27% of women under the age 
of 35 engage in indoor tanning in the United States [18]. Among 
genetic factors associated with melanoma, polymorphisms in the 
melanocortin- 1 receptor (MCIR) are common and exhibit 
moderate penetrance. MCIR is a key regulator of melanin 
synthesis, playing a major role in hair and skin pigmentation [19]. 
MCIR is also involved in UV-induced DNA repair [20]; and 
several of its 80 known variants are associated with melanoma risk 
[19,21]. 

What has yet to be determined is whether genetic markers and 
indoor UV exposure can improve the performance of a clinically- 
based melanoma risk prediction model. Our primary objective of 
this study, therefore, was to determine if adding a genetic marker 
(e.g. MCIR genotype) and indoor UV exposure measures to a 
cUnically-based melanoma risk assessment model could improve its 
predictive ability. A secondary objective was to investigate the 
model's performance when stratified by patient subpopulations 
defined by nevus phenotype, as a way to assess the possible effect 
of melanoma heterogeneity in risk prediction models. 

Materials and Methods 

Study Subjects 

The current study, based on information collected from the 
Minnesota Skin Health Study (SHS), and was deemed exempt 
from review by the IRB at New York University Langone Medical 
Center because a) our study consisted of existing data, and b) these 
data did not include personal identifying information. 

The data analyzed in this study are based on information 
collected from a subset of subjects from the SHS, a study approved 
by the Institutional Re\'iew Board (IRB) at the University of 
Minnesota. This population-based case control study evaluated the 
association between outdoor and indoor UV exposure and 
melanoma risk [11]. Briefly, cases were accrued through the 
Minnesota Cancer Surveillance System, the state's cancer registry. 
Indi\iduals' aged 25 to 59 diagnosed with in\asive cutaneous 
melanoma, of any histologic subtype, betwcx-n July 2004 and 
December 2007 were eligible for enrollment. Controls were 
randomly selected from the Minnesota state driver's license list 
and matched to cases in a 1 : 1 ratio on age (in 5-year age groups) 
and gender. Cases and controls were required to be English- 
speaking and to have a telephone number. In addition to 
providing mouthwash samples, subjects received a self-adminis- 
tered questionnaire, from which selected information was used to 
facilitate a detailed 1-hour telephone interview that collected 
information pertaining to indoor tanning, outdoor sun exposure 
measures, sunscreen use, family history of melanoma, and host 
characteristics (hair, eye, and skin color; freckling; and mole 
phenotype). Nevus/mole phenotype was assessed using cartoon 
diagrams illustrating 4 categories of nevus density ("none", "few". 



"some", or "many" moles). For the purpose of this study, subjects 
who reported having "none" or "few" moles were categorized as 
"nevus-resistant", while those with "some" or "many" moles were 
categorized as "nevus-prone" (Figure SI). Further description of 
materials and methods from the SHS, including assessment of bias 
and questionnaire development, can be found in the original 
publication [1 1]. 

MCIR Genotyping 

Based on availability of germ line DNA, extracted from 
mouthwash .specimens, our analytic sample was comprised of 
1640 subjects from the SHS. MCIR genotyping was performed 
via Sanger sequencing of the entire coding region of MCIR. 
Variants were detected using Mutation Surveyor Software and 
confirmed by visual inspection. 

We based the categorization of MCIR variants on the results of 
a meta-analysis by Whiteman et al. [13]. Variants strongly 
associated with red hair color (RHC) included D84E, R151C, 
R160W, and D294H were designated as "R". Variants weakly 
associated with red hair color (NRHC) included V60L, V92M, 
R142H, I155T, and R163Q,were denoted "r". [13,14]. Histor- 
ically, studies have not addressed the appropriate assignment of 
rare variants and insertions/deletions (indels) as "R" or "r" 
[13,14,22]. We assigned indels and rare variants as "R" or "r" 
based on their association with melanoma case-control status using 
a two degree-of-freedom chi-square test (Pearson or Fisher Exact) 
for genotypes. We also assessed their melanoma odds ratios when 
combined into genotypes with consensus, "R", and "r" alleles. 
The specific assignments of these genotypes are described in the 
Results. 

Statistics/Model Building 

To identify the host and environmental exposure characteristics 
collected by the SHS that were most appropriate for inclusion in a 
risk prediction model, we first applied univariate logistic regression 
model on each measure to obtain the unadjusted odds ratio and 
the corresponding p-value from the Wald test, in order to select 
the characteristics significantly associated with melanoma case- 
control status (p-value<0.05). As the categorical measures of UV 
exposure are all ordinal variables, we applied the same univariate 
logistic regression model regarding each ordinal measure as a 
continuous variable. We also calculated the melanoma odds ratios 
(OR) for each characteristic. Of the measures significandy 
associated with melanoma, those chosen for the development of 
our risk prediction model included a) host characteristics used in 
previously developed melanoma risk models, b) outdoor UV 
exposure measures, and c) indoor UV exposure measures. 

We used multivariate logistic regression to calculate odds ratios 
(OR) and 95% confidence intervals (CI) to determine the degree of 
melanoma risk associated with 1 ) selected host and environmental 
risk factors, 2) and MCIR genotypes. 

The first model (Model A) was based only on host risk factors. 
We created additional models (Models B-E) that included difierent 
combinations of UV exposure measures (outdoor and indoor) 
and/ or MCIR genotype to assess their effect on the predictive 
ability of the clinical model. For Model F, we combined UV 
exposure measures and MC 1 R genotype categories together with 
the host risk factors in Model A. 

We compared the discriminative power of the risk indices using 
the area under the Receiver Operating Characteristic (ROC) 
curve metric. DeLong's test was used to test the significance of the 
incremental increase of AUCs under the ROC curves between 
Models B through F and the baseline Model A. We evaluated the 
models' calibration using the Hosmer-Lemeshow test. The model 
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was cross- validated to obtain the 95 % confidence interval of the 
area under the ROC curve (AUG). Seventy five percent of data 
was randomly chosen as the training sample to develop the model. 
This model, with the same regression coefficients, was then applied 
to the remaining 25% of data to assess the AUG of the model. This 
procedure was repeated 10,000 times, resulting in a sample of 
10,000 estimated AUG's. The 95% confidence interval (GI) of the 
AUG was based on the 2.5% percentile and 97.5% percentile of 
the set of estimated AUGs. We then analyzed model performance 
in subgroups stratified by mole phenotype to assess variation in the 
models' discriminative ability in subjects with and without an 
increased number of nevi. 

Results 

Patient Characteristics and their Association with 
IVlelanoma Risl< 

Phenotypic and Environmental Risk Factors Associated 
with Melanoma. Host information and MGIR sequencing 
data were available for 1640 subjects. Cases and controls had a 
similar mean age and gender distribution. Several host character- 
istics were significantly associated with melanoma risk such as light 
pigmentation, freckhng, and an increased number of nevi (detailed 
results shown in Table 1). Our study group is representative of the 
complete sample of cases and controls from the SHS study, as 
demonstrated in our comparison of participant characteristics in 
Table 1. 

Environmental Risk Factors. We separately analyzed 52 
measures of lifetime and decade-specific sun exposure [23], indoor 
tanning exposure and sunburn history, many of which were 
significandy associated with melanoma risk (Tables SI & S2 in File 
SI) [11]. The incorporation of too many predictor variables could 
over fit and negatively impact the integrity of the model [24] . We 
therefore selected environmental characteristics with the most 
appropriate association with melanoma case-control status, which 
we based on statistical significance and melanoma odds ratios 
(OR). The two factors that met our inclusion criteria were 
"lifetime number of outdoor sunburns", and "frequency (hours) of 
indoor tanning". Subjects with no indoor UV exposure were 
assigned to the category of 0 hours of indoor tanning. 

Indoor tanning frequency was measured several ways (Table S2 
in File SI), nearly all of which were significant on univariate 
analysis. We considered hours of indoor tanning to be the most 
appropriate measure, as it minimizes interpersonal variation 
associated with the number and duration of individual indoor 
tanning sessions. There was a significant dose-response relation- 
ship between hours of indoor tanning and melanoma risk, with 
greater than 10 hours of indoor tanning associated with an 
adjusted OR of 2.32 (95% CI 1.73-3.11; p<0.01) in multivariate 
analysis (Table 1). In contrast, lifetime number of sunburns was 
only associated with a significant increase in melanoma risk in 
univariate analysis, but not multivariate analysis. Greater than 5 
sunburns was associated with an OR of 3.23 (95% CI = 1.81-5.76; 
p<0.01) in univariate analysis and 1.60 (95% CI = 0.85-3.01; 
p = 0.15) multivariate analysis. 

Assignment of MCIR Genotype 

Table 2 displays the variant frequencies of the 9 most common 
MGIR polymorphisms. The remaining non-synonymous poly- 
morphisms were categorized as rare variants. Indels comprised 
their own category. We investigated the most appropriate 
categorization of indels and rare variants by combining them into 
genotypes with consensus, "R" or "r" variants (Table 3). Rare 
variants, when combined into genotypes with "R" or "r" alleles. 
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Table 2. MC1R Variant Allele Frequency* 



MCI R Variants 


Cases N (%) 


Controls N (%) 


P-value** 


Allele Dominance (R/r) 


V60L 


192 (21.94) 


146 (19.08) 


0.17 


'' 


V92M 


138 (15.77) 


129 (16.86) 


0.60 


'' 


R163Q 


98 (11.20) 


94 (12.29) 


0.54 


'' 


R142H 


11 (1.26) 


10 (1.31) 


0.90 




I15ST 


18 (2.06) 


6 (0.78) 


0.05 




D84E 


30 (3.43) 


18 (2.35) 


0.25 


R 


R151C 


205(23.43) 


114 (14.90) 


<0.01 


R 


Rieow 


173 (19.77) 


121 (15.82) 


0.04 


R 


D294H 


46 (5.26) 


19 (2.48) 


0.01 


R 


Indels 


22 (2.51) 


6 (0.78) 


0.01 


R 


Rare 


26 (2.97) 


14 (1.83) 


0.18 


R 



*Allele frequency is determined from the total number of chromosomes genotyped. 
** Chi^ or Fischer's exact test, as appropriate. 
doi:1 0.1 371 /journal.pone.Ol 01 507.t002 



were associated with a melanoma OR of 6.63 (95% CI= 1.89- 
23.27, p<0.01). Due to their strong association with melanoma 
case-control status, we classified rare variants as "R" alleles. 
Similarly, we also classified indels as "R" alleles, as genotypes that 
combined indels with consensus alleles had a melanoma OR of 
3.65 (95% CI = 1 . 14-1 1 .65, p = 0.03), and indels in genotypes with 
"R", "r", or "rare" had an OR of 7.29 (95% CI = 1.59-33.34, 
p = 0.01). In Table 4, we show the melanoma odds ratios of the 
assigned MCIR genotypes. Subjects with "r/r" or genotypes 
containing "R" alleles had statistically significantly higher risk of 
melanoma compared to subjects with "consensus" genotypes, with 
the "R/R" genotype carrying the greatest melanoma risk with an 
OR of 4.31 (95% CI; 2.69-6.89; p<0.01). 

Melanoma Risk Model 

To determine the utility of including MCIR genotype and 
environmental exposure measures (outdoor and indoor UV) in 
melanoma risk prediction, we first developed a model based on 
well-established patient phenotypic factors and age (Table 5, 
Model A), which produced an AUG of 0.72 (95% CI 0.69-0.74). 
We then separately added outdoor and indoor UV exposure to the 
model to assess their individual contribution to its predictive ability 
(Models B and C, respectively). As seen in Table 5, outdoor UV 
exposure did not significandy improve the model's performance 
(AUC 0.72, 95% CI = 0.69-0.74; p = 0.66), while the addition of 
indoor UV exposure resulted in a statistically significant improve- 
ment with an AUC of 0.73 (95% CI = 0.70-0.75; p = 0.03). 



Interestingly, Model D, which includes both outdoor and indoor 
UV exposure, had the same AUC as Model C (AUC = 0.73, 95% 
CI = 0.71-0.75; p = 0.02), suggesting that outdoor UV exposure 
does not contribute to the performance of the model. The addition 
of MC IR genotype to Model A (Model E) resulted in a statistically 
significant improvement compared to the baseline model, with an 
AUC of 0.73 (95% CI = 0.70-0.75; p = 0.02). We interpret this 
improvement in model performance to be a demonstration of 
MClR's functional role outside of hair color (i.e. DNA repair). A 
full model (Model F) combining phenotypic factors, outdoor and 
indoor UV exposure, and MCIR genotype performed the best 
with an AUC of 0.74 (95% CI = 0.71-0.76; p<0.01). These results 
were consistent in cross-validation analysis, and the models were 
determined to be well calibrated, as demonstrated by the Hosmer- 
Lemeshow Goodness-of-fit test with p-values greater than 0.05 
(Table 5). 

Testing the Risk Model in Patient Subsets 

Since one of the most important phenotypic risk factors for 
melanoma is number of nevi, we explored the performance of 
these models in the nevus-resistant and nevus-prone patient 
subsets, as described in the Methods section (Table 6). The AUC 
of Model A (baseline model) decreased in both subsets compared 
to all subjects. The decrease was much more profound in the 
nevus-prone group compared to the nevus-resistant group 
(AUC = 0.60 vs. 0.69). In both groups. Model F (final model) 
performed better than Model A. This incremental improvement 



Table 3. Melanoma Risk of Genotypes Combining Indels and Rare Variants with Conventional IV1C1R Variants. 





Genotype 


Cases (n = 875) 


Controls (n = 765) 


Unadjusted OR (95% CI) 


P-value 


Consensus/consensus 


163 


216 


Ref 




rare/consensus 


10 


11 


1.21 (0.50-2.91) 


0.68 


Indel/consensus 


11 


4 


3.65 (1.14-11.65) 


0.03 


rare/(R or r) 


15 


3 


6.63 (1.89-23.27) 


<0.01 


lndel/(r, R, or rare) 


11 


2 


7.29 (1.59-33.34) 


0.01 



"R": D84E, R151C, R160W, D294H. 
"r": V60L, V92M, R163Q, R142H, I155T. 
doi:l 0.1 371/journal.pone.Ol 01 507.t003 
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Table 4. Melanoma Risk Based 


on MCIR Genotyping including Indels 


and Rare Variants. 






Genotype 


Cases (n = 875) 


Controls (n = 765) 


Unadjusted OR (95% CI) 


P-value 


Consensus/consensus 


163 


216 




Ref 




r/consensus 


203 


220 




1.22 (0.93-1.62) 


0.18 


r/r 


67 


55 




1.61 (1.07-2.43) 


0.03 


R/Consensus 


208 


172 




1.60 (1.20-2.13) 


<0.01 


R/r 


143 


74 




2.56 (1.81-3.62) 


<0.01 


R/R 


91 


28 




4.31 (2.69-6.89) 


<0.01 



"R": D84E, R151C, R160W, D294H, indels, rare variants, 
"r": V60L, V92M, R163Q, R142H, I155T. 
doi:1 0.1 371 /journal.pone.Ol 01 507.t004 



was somewhat larger in the nevus-prone group with an increase in 
AUG from 0.60 to 0.67, compared to the nevus-resistant subgroup 
with an increase in AUG from 0.69 to 0.72. In analyzing the step- 
wise improvement of adding UV exposure measures to the 
baseline Model A, the role of indoor UV exposure (Model C) 
appears to have a more profound impact on model performance in 
the nevus-prone subgroup with an increase in AUG from 0.60 to 
0.65 compared to the nevus-resistant subgroup with an increase in 
AUG from 0.69 to 0.71. 

Discussion 

Our study has several key findings. First, we found that indels 
and rare variants were associated with a significantly increased 
melanoma risk. Second, we demonstrated that adding MCIR 
genotype and indoor UV exposure data to a phenotypic 
melanoma risk model results in a small, but statistically significant 
increase in predictive ability, which was upheld when hair color 
measures were removed from the model (data not shown). The 
incremental increase in AUG corresponding to adding the 
predictors is quantitatively small, but it is well known that 
incremental AUG is a conservative measure of discrimination 
improvement [25-27]. These findings support the potential utility 
of genetic risk markers to improve the recognition of the more 



Table 5. Description of Melanoma Risk Models. 



than 50% of melanoma patients that lack common phenotypic risk 
factors [5,6,1 1,28]. Lastly, we observed substantial variation in the 
contribution of indoor UV exposure to the model's performance 
when subjects were stratified by mole phenotype. This stratifica- 
tion is further supported by interaction analyses suggesting that 
indoor UV exposure may confer differing increases in melanoma 
risk between these subgroups (p = 0.06, unpublished data). This 
suggests that the inherited genetic variants that contribute to the 
"nevus prone" and "nevus resistant" phenotypes may interact 
differently with UV exposure to affect melanoma risk. 

Whiteman et al. first described the concept that different mole 
phenotypes may be associated with different melanoma causal 
pathways as part of "The Divergent Pathway Theory" [29]. This 
theory provides a conceptual framework connecting the epidemi- 
ologic heterogeneity of melanoma with melanoma risk factors, 
particularly mole phenotype [30]. Published reports demonstrate 
that nevogenesis has a strong genetic component [9,10]; therefore, 
nevus-prone and nevus-resistant individuals likely have different 
germline genetic variants that influence their mole phenotypes and 
melanoma risk. Since the number of nevi is extremely important in 
melanoma risk prediction, stratification by nevus phenotype 
allowed us to assess whether the additional risk factors (indoor 
and outdoor UV measures, and MGIR genotype) could improve 



Hosmer-Lemeshow 



Model 


Variables 


All Subjects (N 


= 1640) 


Cross Validation 


Goodness-of-Fit test 






AUC (95% CI) 


P value* 


AUC (95% CI) 


P value 


A 


Age 


0.72 (0.69-0.74) 


Ref 


0.70 (0.66-0.75) 


0.90 


Gender 


Hair color 


Eye color 


Skin color 


Freckling 


Mole phenotype 


B 


A+total sunburns 


0.72 (0.69-0.74) 


0.66 


0.70 (0.66-0.75) 


0.60 


C 


A+hours of indoor tanning 


0.73 (0.70-0.75) 


0.03 


0.71 (0.67-0.76) 


0.25 


D 


A+total sunburns+hours of indoor tanning 0.73 (0.71-0.75) 


0.02 


0.71 (0.67-0.75) 


0.08 


E 


A+MC1 R Genotype 


0.73 (0.70-0.75) 


0.02 


0.71 (0.66-0.75) 


0.69 


F (Full model) 


A+D+E 


0.74 (0.71-0.76) 


<0.01 


0.72 (0.67-0.76) 


0.65 



* p-values in reference to Model A. 
doi:l 0.1 371/journal.pone.Ol 01 507.t005 
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Table 6. Melanoma Risk Models Stratified by Mole Phenotype. 



Model 


All patients (N = 1640) 




Nevus Resistant* 


(N = 1321) 


Nevus Prone* (N 


= 319) 




AUC (95% CI) 


P value 


AUC (95% CI) 


P value 


AUC (95% CI) 


P value 


A 


0.72 (0.69-0.74) 


Ref 


0.69 (0.67-0.72) 


Ref 


0.60 (0.53-0.68) 


Ref 


B 


0.72 (0.69-0.74) 


0.66 


0.70 (0.67-0.72) 


0.50 


0.61 (0.53-0.68) 


0.44 


C 


0.73 (0.70-0.75) 


0.03 


0.71 (0.68-0.73) 


0.09 


0.65 (0.57-0.72) 


0.03 


D 


0.73 (0.71-0.75) 


0.02 


0.71 (0.68-0.74) 


0.06 


0.65 (0.58-0.72) 


0.02 


E 


0.73 (0.70-0.75) 


0.014 


0.70 (0.68-0.73) 


0.03 


0.62 (0.55-0.70) 


0.15 


F (Full model) 


0.74 (0.71-0.76) 


<0.01 


0.72 (0.69-0.74) 


0.01 


0.67 (0.60-0.74) 


0.01 



*The AUC for Nevus resistant and Nevus prone is calculated based on the coefficients obtained from tlie "All Patient" model. 
doi:l 0.1 371 /journal.pone.Ol 01 507.t006 



the model's performance among the lower-risk, nevus-resistant 
individuals who comprised 72% of the melanoma cases. Although 
the final model (Model F) performed better in the nevus-resistant 
compared to the nevus-prone subgroups (AUC = 0.72 vs. 0.67), 
the degree of improvement from Model A to Model F was 
somewhat larger in the nevus-prone group. While we recognize 
die inherent limitations of subgroup analyses, this variation in 
model performance with the addition of indoor UV exposure 
provides support for further studies investigating the interaction 
between environmental exposure and an individual's genetic 
propensity to develop nevi [31]. 

Growing evidence supports the utility of targeted cancer 
screening. In lung cancer, screening a high-risk population led to 
a significant decrease in mortality [4,32]. Similarly, a population- 
based melanoma screening program in Germany was associated 
with a nearly 50% reduction in melanoma mortality [3]; however, 
we estimate the number needed to screen to prevent one 
melanoma death at approximately 127,000. This potentially high 
cost of population-based melanoma screening could be reduced 
through a targeted approach similar to that used in lung cancer. 
The risk assessment model we describe offers one approach to 
developing a tool to identify individuals for targeted screening, 
particularly the potential benefit of using genetic information. 
Adding MCIR genotype resulted in a small, yet significant 
improvement in the predictive ability of our model. It is also worth 
noting that this improvement was based on a single gene. In 
comparison, risk indices for breast and prostate cancers require 
several genetic markers to produce increases of similar magnitude 
[22,24,28,33-36]. 

MCIR genotype has been incorporated into two other 
preliminary risk models [37-39]. A Greek hospital-based study 
incorporated 8 melanoma-related single nucleotide polymor- 
phisms (SNP's), some of which were MCIR variants, into a 
clinically-based risk model, but investigators did not find a 
significant improvement in their clinical model's performance 
with the addition of the SNP's [38] . Most recently. Gust et al. 
demonstrated that adding MC 1 R genotypic information increased 
in the performance of a baseline melanoma risk prediction model, 
supporting our findings [39] . Of note, their baseline model used 
age, sex, city of recruitment, and self-reported European ancestry 
as covariates, which are very different from our baseline model, 
suggesting that MGIR genotype may be a robust factor to help 
identify patients at increased risk for melanoma. 

Our preliminary melanoma risk model has several strengths 
compared to existing models [40-43]. First, we developed the 
model using a high-quality, population-based case-control study of 
over 1600 subjects. Second, we noted a small, yet significant 



increase in the performance of our model with the addition of UV 
exposure data, particularly indoor tanning. Third, we character- 
ized each patient for their MCIR genotype, rather than the 
presence of specific variants. We believe this is a more 
comprehensive approach as it accounts for both alleles. Finally, 
we are the first to demonstrate a variation in model performance 
when stratified by mole phenotype, which suggests that more than 
one melanoma risk model may be needed to address melanoma 
etiologic heterogeneity. 

In addition to our model's strengths, there are also limitations. 
The development of a successful melanoma risk model using self- 
reported host and environmental exposure factors has the benefit 
of being easily accessible to the general population and potentially 
cost-effective; conversely, the potential for inaccurate self-assess- 
ments of host characteristics (i.e. number of nevi), and recall bias 
with respect to UV exposure may result in misclassification of 
melanoma risk. Of note, the risk of recall bias associated with self- 
reported risk factors was addressed in the parent study. Such bias 
was analyzed and found not to influence the odds ratios for various 
phenotypic factors [11]. Secondly, our study (and the parent study) 
excluded subjects older than age 59 due to the decreasing 
prevalence of indoor tanning among older patient cohorts. 
Subsequent studies are needed to test the model in populations 
that include older individuals. Finally, the small sample size of 
nevus-prone subjects (n = 319), is a potential limitation of our 
subgroup analysis, and will need to be validated in a larger sample 
size. The performance of modified MCIR genotyping may also 
vary by ancestry and/ or geographic locations. 

Understandably, our melanoma risk model (Model F) requires 
replication in additional patient cohorts; however, as a hypothesis- 
generating model, it is promising that variations in a single gene 
can significantly improve the predictive performance of a model 
incorporating host and environmental measures. Very recently. 
Fang et al. demonstrated that the addition of 1 1 SNP's identified 
in melanoma GWAS studies to a basic phenotypic model (i.e. age, 
gender, hair color, eye color, and skin color) resulted in a 7.8% 
increase in AUC from 0.64 to 0.69 [44]. This fmding supports the 
role of genetic markers to potentially improve melanoma risk 
prediction, and suggests that adding these and/or additional 
genetic markers may further improve the performance of our risk 
model. 

In conclusion, we demonstrated that the inclusion of indoor 
tanning measures and MCIR genotype improve the predictive 
ability of a clinically-based melanoma risk model. Variation in the 
risk model's discriminative ability when applied to phenotypic 
subgroups suggests that the influence of certain melanoma risk 
factors may vary by a patient's clinical characteristics, supporting 
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the disease heterogeneity of melanoma as defined by the Divergent 
Pathway Theory. Finally, the performance improvement by 
including MCIR genotypes supports existing evidence that 
genetics can be used to improve melanoma risk prediction. 

Supporting Information 

Figure SI Diagrams illustrating 4 categories of self- 
reported nevus density. Clockwise from top, diagrams 1-4 
correspond with having "none", "few", "some", or "many" 
moles, respectively. 
(TIF) 

References 

1. Rics L, Mclbcrt D, Krapcho M, Stinchcomb DG, Howladcr N, ct al. (2008) 
SEER Cancer Statistics Review, 1975-2005. 

2. Gimotty PA, Elder DE, Fraker DL, Botbyl J, Sellers K, ct al. (2007) 
Identification of high-risk patients among those diagnosed with thin cutaneous 
melanomas. J Clin Oncol 25: 1129-1134. 

3. Katalinic A, Waldmann A, Weinstock MA, Geller AC, Eisemann N, et al. (20 1 2) 
Does skin cancer screening save lives?: An observational study comparing trends 
in melanoma mortality in regions with and without screening. Cancer. 

4. Scstini P (201 1) Reduced lung-cancer mortality with CT screening. N Engl J Med 
,365: 2037; author reply 2037-20,38. 

5. Tucker MA, Halpcrn A, Holly EA, Hartge P, Elder DE, ct al. (1997) Clinically 
recognized dysplastic nevi. A central risk factor for cutaneous melanoma. JAMA 
277: 1439-1444. 

6. Gandini S, Sera F, Cattaruzza MS, Pasquini P, Abeni D, et al. (2005) Meta- 
analysis of risk factors for cutaneous melanoma: I. Common and atypical naevi. 
Eur J Cancer 41: 28^4. 

7. Gandini S, Sera F, Cattaruzza MS, Pasquini P, Zanetti R, et al. (2005) Meta- 
analysis of risk factors for cutaneous melanoma: 111. Family history, actinic 
damage and phenot^'pic factors. Eur J Cancer 41; 2040-2059. 

8. Whitemaii DC, Pavan VVJ, Bastian BC (2011) The melanomas; a synthesis of 
epidemiological, clinical, histopathological, genetic, and biological aspects, 
supporting distinct subtypes, causal pathways, and cells of origin. Pigment Cell 
Melanoma Res 24: 879-897. 

9. Newton-Bishop JA, Chang YM, Bes MM, Taylor JC, Bakker B, et al. (2010) 
Melanocytic nevi, nevus genes, and melanoma risk in a large case-control study 
in the United Kingdom. Cancer Epidemiol Biomarkers Prev 19: 2043—2054. 

10. Falchi M, Bataille V, Hayward NK, Dull> DL. Bishop JA, et al. (2009) Genome- 
wide association study identifies variants at 9p21 and 22ql3 associated with 
development of cutaneous nevi. Nat (ienet 41; 915-919. 

11. Lazovich D, Vogcl RI, Berwick M, Weinstock MA, Anderson KE, et al. (2010) 
Indoor tanning and risk of melanoma: a case-control study in a h^hly exposed 
population. Cancer Epidemiol Biomarkers Prev 19: 1557—1568. 

12. Iharrola-Villava M, Hu HH, Guedj M, Fernandez LP, Descamps V, et al. (2012) 
MCIR, SLC45.^2 and LYR genetic variants involved in melanoma suscepti- 
bility in Southern Luropean populations; Results from a Meta-analysis. 
Eur J Cancer. 

13. WiUiams PF, Olscn CM, Hayward NK, Whiteman DC (2011) Melanocortin 1 
receptor and risk of cutaneous melanoma; a meta-analysis and estimates of 
population burden. IntJ Cancer 129: 1730-1740. 

14. Kanetsky PA, Panossian S, Elder DE, Guerry D, Ming ME, et al. (2010) Does 
MCIR genotype convey information about melanoma risk beyond risk 
phenotypes? Cancer 116: 2416-2428. 

15. Hill VK, Gartner JJ, Samuels Y, Goldstein A.\I (2013) The genetics of 
melanoma; recent advances. Annu Rev (Genomics Hum Genet 14; 257—279. 

16. Boniol M, Autier P, Boyle P, (iandiiii S (2012) (?.utaneous melanoma 
attributable to sunbed use; systematic review and meta-analysis. BMJ 345; 
e4757. 

17. M B, P A, P B, S G (2012) Cutaneous melanoma attributable to sunbed use: 
systematic review and meta-analysis. 345: e8503. 

18. Guy GP, Berkowitz Z, Watson M, Holman DM, Richardson LC (2013) Indoor 
tanning among young non-Hispanic white females. JAMA Intern Med 173; 
1920-1922. 

19. Fargnoli MC, Gandini S, Peris K, Maisonneuve P, Raimondi S (2010) MCIR 
variants increase melanoma risk in families with CDKN2A mutations: a meta- 
analysis. Eur J Cancer 46; 1413-1420. 

20. Wong SS, Ainger SA, LeonardJH, Sturm RA (2012) MCIR variant allele effects 
on UVR-induced phosphorylation of p38, p53, and DDB2 repair protein 
responses in melanocytic cells in culture. J Invest Dermatol 132: 1452—1461. 

21. Raimondi S, Sera F, Gandini S, lodice S, Caini S, et al. (2008) MCIR variants, 
melanoma and red hair color phcnot^'pe; a meta-analysis. IntJ Cancer 122: 
275,3-2760. 

22. Elfakir A, Ezzedinc K, Latreillc J, Ambroisine L, Jdid R, et al. (2010) Functional 
MClR-gene variants are associated with increased risk for severe photoaging of 
facial skin. J Invest Dermatol 130; 1 107-1 1 15. 



File SI Contains Tables SI and S2. Table SI. Univariate 
Analysis of Association between Melanoma and Outdoor UV 
Exposure Variables Collected by the SHS. Table S2. Univariate 

Analysis of Association between Melanoma and the (categorical) 

Indoor tanning Variables Collected by the SHS. 

(DOCX) 

Author Contributions 

Conccix'cd and designed the experiments: MB DL DP LP YS. Performed 
the experiments: MQ_EZ EN LP. Analyzed the data: MQ,EZ YS LP DP. 
Contributed rcagents/materials/analysis tools: MB DL. Wrote the paper: 
LP MQ,YS MB DL DP. 



23. Kricker A, Vajdic CM, Armstrong BK (2005) Reliability and validity of a 
telephone c|uestioniiaire for estimating lifetime personal sun exposure in 
epidemiologic studies. Cancer Epidemiol Biomarkers Prev 14: 2427—2432. 

24. Ranucci M, Castelvecchio S, Menicanti L, Frigiola A, Pelissero G (2010) 
Accuracy, calibration and clinical performance of the EuroSCORE: can we 
reduce the number of variables? Eur J Cardiothorac Surg 37: 724—729. 

25. Vickers AJ, Cronin AM, Begg CB (2011) One statistical test is sufficient for 
assessing new predictive markers. BMC Med Res Methodol 11: 13. 

26. Gail MH (2008) Discriminatory aceuracv from single-nueleotidc polymorphisms 
in models to predict breast cancer risk. J Natl Cancer Inst 100; 1037—1041. 

27. Folsom AR, Chambless LE, Ballantvnc CM, CorcshJ, Heiss G, et al. (2006) An 
assessment of incremental coronary risk prediction using C-reactive protein and 
other novel risk markers; the atherosclerosis risk in communities study. Arch 
Intern Med 166; 1368-1373. 

28. Udayakumar D, Tsao H (2009) Moderate- to low-risk variant alleles of 
cutaneous mahgnancies and nevi: lessons from genome-wide association studies. 
Genome Med I: 95. 

29. Olsen CM, Zens MS, Sttikel TA, Sacerdote C, Chang YM, et al. (2009) Nevus 
density and melanoma risk in women: a pooled analysis to test the divergent 
pathway hj'jiothcsis. IntJ Cancer 124: 937-944. 

30. Whiteman DC (2010) Testing the divergent pathway hypothesis for melanoma; 
recent findings and future challenges. Expert Rev Anticancer 'I her 10; 615-618. 

31. Whiteman DC, Brown RM, Purdie DM, Hughes MG (2005) Melanocytic nevi 
in very young children; the role of phenotype, sun exposure, and sun protection. 
J Am Acad Dermatol 52: 40-47. 

32. Smith RA, Cokkinides V, Brawley OW (2009) Cancer screening in the United 
States, 2009: a review of current American Cancer Society guidelines and issues 
in cancer screening. CA Cancer J Clin 59: 27—41. 

33. Mealiffc ME, Stokowski RP, Rhccs BK, Prentice RL, Pettinger M, ct al. (2010) 
Assessment of clinical validity of a breast cancer risk model combining genetic 
and clinical information. J Natl Cancer Inst 102; 1618—1627. 

34. Lindstrom S, Schumacher FR, Cox D, Travis RC, ."Ubanes D, ct al. (2012) 
Common genetic variants in prostate cancer risk prediction-results from the 
NCI Breast and Prostate Cancer Cohort Consortium (BPC3). Cancer Epidemiol 
Biomarkers Prev 21; 437-444. 

35. Macinnis RJ, Antoniou AC, Eeles RA, Severi G, Al Olama AA, et al. (201 1) A 
risk prediction algorithm based on family history and common genetic variants: 
application to prostate cancer with potential clinical impact. Genet Epidemiol 
35; ,549-5,56. 

36. Chatterjee N, ParkJH, Caporaso N, Ciail MH (2011) Predicting the future of 
genetic risk prediction. Cancer Epidemiol Biomarkers Prev 20; 3—8. 

37. Whiteman DC, (ireen AC (200.5) A risk prediction tool for melanoma? Cancer 
Epidemiol Biomarkers Prev 14: 761—763. 

38. Stefanaki I, Panagiotou OA, Kodela E, Gogas H, Kypreou KP, et al. (2013) 
Rephcation and predictive value of SNPs associated with melanoma and 
pigmentation traits in a Southern European case-control study. PLoS One 8: 
e55712. 

39. Gust AE, Goumas C, Vuong K, Davies JR, Barrett JH, et al. (2013) MCIR 
genotype as a predictor of early-onset melanoma, compared with self-reported 
and physician-measured traditional risk factors; an Australian casc-control- 
family study. BMC Cancer 13; 406. 

40. Fortes C, Mastroeni S, Bakos L, Antonelli G, Alessandroni L, et al. (2010) 
Identifying individuals at h^h risk of melanoma: a simple tool. Eur J Cancer 
Prev 19: 393-400. 

41. Cho E, Rosner BA, Feskanich D, Colditz GA (2005) Risk factors and individual 
probabilities of melanoma for whites. J Clin Oncol 23: 2669-2675. 

42. Fears TR, Guerry Dt, Pfeiffer RM, Sagebiel RW, Elder DE, et al. (2006) 
Identifying individuals at high risk of melanoma: a practical predictor of absolute 
risk. J Clin Oncol 24; 3590-3596. 

43. Mar V, Wolfe R, Kelly JW (201 1) Predicting melanoma risk for the Austrahan 
population. AustralasJ Dermatol 52; 109—116. 

44. Fang S. ILuiJ, Zhang .\L Wang LL. Wei Q, et al. {2013)Joint Effect of Multiple 
Common SNPs Predicts Melanoma Susceptibility. PLoS One 8; e85642. 



PLOS ONE I www.plosone.org 



8 



July 2014 I Volume 9 | Issue 7 | e101507 



